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Abstract 

The internal motions of DNA immersed in bio-fluid are investigated. The 
interactions between the fragments of DNA and the surrounding bio-fluid are 
modeled using the gauge fluid lagrangian. In the model, the bio-fluid is cou- 
pled to the standard gauge invariant bosonic lagrangian describing the DNA. 
It is shown that at non-relativistic limit various equation of motions, from the 
well-known Sine-Gordon equation to the simultaneous nonlinear equations, 
can be constructed within a single framework. The effects of bio-fluid are in- 
vestigated for two cases : single and double stranded DNA. It is argued that 
the small and large amplitudes of a single stranded DNA motion immersed 
in bio-fluid can be explained in a natural way within the model as a solitonic 
wave regardless with the fluid velocity. In contrary the double stranded DNA 
behaves as regular or damped harmonic oscillator and is highly depending on 
the fluid velocity. 
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1 Introduction 



Both deoxyribo- and ribo-nucleic acid (DNA and RNA) have been recognized as 
the most important biomolecules. Especially DNA helical structures undergo a 

very complex dynamics which plays several important roles in various biological 
phenomena such as storage of information, inheritance (replication, etc) and the 
usage of genetic information (transcription, etc). The importance of biopolymers 
hke DNA/RNA is motivated by estabhshed observations that the homologous re- 
combination is preceded by recognition and local pairing of intact double stranded 
DNA fragments, rather than involving known recombination proteins. Therefore, it 
should be attributed to direct DNA-DNA interactions whose physical origin has not 
been understood [1, 2]. Experimentally, the physical properties of DNA/RNA have 
been measured in many works, for example : the DNA single-molecule [3, 4, 5], dou- 
ble stranded DNA forming bubbles [6], the DNA/RNA nucleoside and nucleotides 
[7], the structural transitions of DNA through torques measurements [8], the ther- 
modynamic fluctuations of DNA in a reacting system [9], the stretching DNA with 
a receding meniscus [10], the electrical transport through single DNA molecules [11] 
and so forth. 

Prom physical point of view, a biopolymer like DNA molecule is considered as 
a system consisting of many interacting matters in a particular configuration of 
space-time. Some models treats this kind of DNA dynamics as the phenomena 
of nonlinear excitations like soliton. This type of models has been pioneered by 
Englander et.al. using nonlinear dynamics relevant to the transcription process in 
terms of coupled pendulum chain which generates the sine-Gordon equation and 
its classical solitons [12]. Further, Davydov described the alpha helices in quantum 
solitons [13]. Following these suggestions, a number of models for the nonlinear DNA 
have been elaborated in the last decades, in both classical or quantum approaches 
[14, 15, 16]. A typical classical approach is the so-called PDB model which takes into 
account twisted DNA molecules [17, 18, 19]. On the other hand, there are several 
models based on the particle interactions [20, 21, 22]. Also, the polyelectrolyte 
model which treats DNA molecule as a cylinder with a net charge homogeneously 
distributed along its surface, and has further been modified to be the electrostatic 
zipper motif for DNA aggregation [23], to solve high dependency of the electrostatic 
interaction between DNA duplexer on surface charge patterns [24]. 

It has also been shown that under particular external conditions the DNA molecules 
form a double helix, and its (transverse, longitudinal and torsional) motions can be 
divided into two main regions : the small and large amplitude of internal motions 
[25] . The small amplitude of motion can be described by the hamiltonian of harmonic 
oscillator. On the other hand, the large amplitude is described by a non-harmonic 
one [26]. Recently, many works have discussed and arrived at the conclusion that 
the large amplitude of internal motion can be considered as a nonlinear dynamical 
system where solitary conformational waves can be excited [14]. Then nonlinear in- 
teraction between molecules in DNA gives rise to a very stable excitation as sohton 
[26, 27]. 

As mentioned above, DNA is not motionless. It is in a constantly wriggling dy- 
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namics state in a medium of bio-organic fluid in the nucleus cell [28]. However, the 
motion of DNA surrounded by fluid is rarely studied. Previous studies are usually 
done by solving the fluid equations and its wave equations simultaneously using 
appropriate boundary conditions. On the other hand, in the Hamiltonian formula- 
tion the viscous force is considered to be comparable with other forces arising from 
Hamiltonian [29, 30]. The solution is then obtained by expansion and performing 
order-by-order calculation. In these approaches, anyway the picture of interaction 
between DNA and its surrounding fluid is not clear. Also, in most models the over- 
damped DNA dynamics are treated by putting some additional terms by hand in 
the differential equation to obtain the non-homogenous ones [31]. The stochastic 
simulations of DNA in flow has been done for a fully parametrized headspring chain 
model by taking into account the fluctuating hydrodynamic interactions [32]. 

In this paper, a new model to describe various internal motions of DNA inspired 
by gauge fluid theory is proposed. The DNA dynamics is modeled as the result 
of interactions among matters in a fluid medium using the relativistic and gauge 
invariant fluid lagrangian. Although the theory is a relativistic one, we take its 
non-relativistic limit at the flnal stage to deal with problems in DNA as done in 
some previous works, for instance in some models using the ideal gas approximation 
[33]. Moreover, the lagrangian is intended for physics at scale of order transport 
mean free paths, that is the transition region where neither a hydrodynamics nor 
kinetic theory is valid. Therefore it fits the current interest of modeling "elementary" 
biomatters like DNA. Just to mention, the lagrangian is originally devoted to model 
the quark gluon plasma (QGP) as a relativistic fluid system [34, 35, 36, 37], inspired 
by the similarity between the dynamical properties of fluid and electromagnetic field 
[38, 39]. The DNA is treated as strongly coupled system like non-Abelian plasmas 
where neither a hydrodynamics nor kinetic theory is really valid. Within the model, 
a single and double stranded biopolymcrs arc described in a general way as the 
results of interactions among the fiuid and matter fields. 

We show in two specific cases how to derive the equation of motion (EOM) and 
investigate the internal motions through its solutions and behaviors. From the EOM 
of a DNA as a single bulk, we argue that small and large amplitude regions of the 
internal motion of DNA are determined by its internal dynamics and interactions 
with surrounding fiuid. On the other hand, in the case of double stranded DNA the 
EOM is solved analytically to investigate the effects of fiuid velocity to its internal 
motion. 

The paper is organized as follows. First we briefiy introduce the theory of gauge 
invariant fiuid lagrangian, and then provide the allowed interactions within the 
model. After explaining how to model DNA using the interactions in the lagrangian, 
we provide two typical examples : 1) the Abelian U(l) case to model the dynamics 
of a single bulk of DNA, and 2) the non-Abelian SU(2) case to describe the internal 
motion of double stranded DNA. Finally, the paper is ended with summary and 
discussion. 
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2 Theoretical background 



Here, a new approach to investigate the interaction between biopolymer and its 
surrounding bio-fluid is discussed using the lagrangian method. Rather putting it 
by hand, the interaction is described in a more natural way from first principle, i.e. 
by introducing some symmetries in the lagrangian under consideration. 



2.1 The lagrangian 

The model is an extension of the original model based on the U(l) gauge theory 
devoted for QGP as a magnctofiuid system [34, 35, 36]. Thereafter it has been 
extended to the non-Abclian case to accommodate a system with many matters, 
either bosonic or fermionic ones [37]. Concerning the fact that an (elementary) 
matter has no intrinsic degree of freedom like spin, it is considerable to represent its 
elementary constituents as scalar (boson) fields governed by the bosonic lagrangian. 



r-m.tter = (5^*)^ (d^^) + y($) , 

where V{^) is the potential. For example in the typical theory, 

y($) = _l^|$t$_i.A($t$)2^ 



(i: 



(2) 



where m$ and A are the mass of matter and the dimensionless coupling constant of 
matter self-interaction. The hermite conjugate is = ($*)'^ for a general complex 
field $. 

We impose the above bosonic lagrangian to be gauge invariant under local 
(in general non-Abehan) gauge transformation [40, 41], [/ = exp[—iT"-9"-{x)] ^ 
1 — iT"-9"-{x) with 9°- <^ 1. T^'s are generators belong to a particular Lie group 
and satisfy certain commutation relation [T",T''] = ijabcj^c ^j^j^ jabc jg ^j^g gjiH- 

symmetric structure constant [42]. The matter field is then transformed as $ 

$' = exp[— iT"^^"(a;)] $, with T" are n x n matrices while $ is an n x 1 multiplet 

containing n elements, i.e. 



$2 



and = ($1 $2 



*n), 



(3) 



for n dimension Lie groups as SU(n), 0{n + 1), etc. It is well-known that the sym- 
metry in Eq. (1) is revealed by introducing gauge fields A^^ which are transformed 

as U^' = [/; - l{df,9'') + f'''"'9''U^, and replacing the derivative with the 

covariant one, = 9^ + igT'^U'j^. Anyway, the number of generators, and also 
gauge bosons, is determined by the dimension of group under consideration. For an 
SU(n) group one has — 1 generators and the index a runs over 1, 2, • • • , — 1. 
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For example the SU(2) group is realized by 2 x 2 matrices T" = |cr" with cr" are the 
Pauli matrices [42], 

1 / 01\ 2 3 f 1 



r — C" QCMi^ (f,\ 



In particular, the Abelian U(l) case is revealed by putting T"'9"-{x) — >■ 9{x), i.e. the 
phase transformation, respectively. 

Finally, the gauge invariancc leads to the total lagrangian with some additional 
terms in the lagrangian to keep its gauge invariance, 

-^matter ~l~ -^gauge ~l~ -^int i (5) 

where, 

1 

4' 

Ant = -gJlW + {^^T''T''^)UIU^'" . (7) 

The strength tensor is -S^j, = d^JJ^ — + gf"'^'^U^^U^, while the 4-vector current 
is, 

j; = -i[{d^^)^n-^^T\d^^)\ . (8) 

The coupling constant g then represents the interaction strength between gauge 
field and matter. We should note that, the current conservation is realized by the 
covariant current J"^" = with J"^" = -i [(P^$)tr"$ - [37]. 

The gauge boson tf^ is interpreted as a "fluid fleld" with velocity u^j,, and takes 
the form [34, 35, 36, 37], 

c/; = (t/«,u«)^«;;0, (9) 

with, 

K^^7'^(l,-v'^), (10) 

where is an auxiliary boson field, while 7" = (1 — |v"p) Here we adopt the 
natural unit, i.e. the light speed c = 1. Eq. (10) is nothing else than rewriting a 
gauge field in terms of its polarization vector and wave function which represents 
the fluid distribution in a system. It has further been shown that the non-relativistic 
fluid equation can be reproduced using Eq. (9) [35, 37]. This fact actually justifles 
us to model the DNA dynamics in a fluid medium using the total lagrangian in Eq. 
(5). 

Now we are ready to model the DNA using the above lagrangian. First, we 
should investigate the allowed interactions in the present theory. 



2.2 The interactions 

In order to be specific, let us consider the potential in Eq. (2) for matter 
lagrangian in Eq. (1). With a complete lagrangian at hand, we can extract m— point 
interactions for fluid and matter with m is the number of relevant legs involved in 
an interaction. We list all allowed interactions below for each element in the matter 
multiplet denoted by the indices i, j. 
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• 2— point interactions : 

The interactions arise through the kinetic and mass terms of matter in Eqs. 
(1) and (2), and the fluid kinetic term in Eq. (6), 

: {d^^*){d^'^i)-^ml^*^i. (11) 
UU : ~ {d^Ut - d^Ul) {dW - d'U"'') . (12) 

• 3— point interactions : 

These interactions are induced by the fluid self-interaction in Eq. (6) and the 
fluid- matter interaction in Eq. (7), 

: i^7^«[(a^$*).$,-$*(5M*)Jt^"''- (13) 
UUU : lgr"'U'^U'''{d^U^-dM^^) , (14) 

• 4— point interactions : 

These interactions are induced through the matter self- interaction in Eq. (1), 
the fluid kinetic term in Eq. (6) and the fluid- matter interaction in Eq. (7), 

$$$$ : -^A , (15) 

UUUU : (16) 
^^UU : g'^^{T''T^)..^jU''^U^'' . (17) 

All of these constitute the so-called Feynman diagrams and its order of mag- 
nitudes that will be used soon in the subsequent section. Now we are ready to 
construct the models relevant for biopolymers. 



3 Modeling the DNA 

Here, we consider two typical examples on how to describe various dynamics of DNA 
within the present model. First we present a model for a single bulk of DNA or a 
fragment of DNA molecule like nucleotide or nucleoside. Further we construct a 
more complicated picture for the double stranded DNA. The model is a new type 
of the mesoscale model of DNA that reduces the complexity of a nucleotide to three 
interactions sites [21]. 

3.1 Single bulk of DNA : the Abelian U(l) model 

The Abelian U(l) lagrangian involves only a single matter and a fluid field. In this 
case, the 3— point interaction in Eq. (14) and the 4— point interaction in Eq. (16) 
vanish. It is also clear that we are not able to construct a reahstic model for a 
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biopolymer composed by several difFerent matters in this case [22]. However, we 
can model the dynamics of a single bulk of DNA or its fragment like nucleoside 
which could be considered as a composite field of sugar and base. This means we 
investigate the internal dynamics of namely DNA molecules through the EOM of 
its fragments and study the basic behaviors. 
The total lagrangian in this case becomes, 

(d^U, - d^U^) {dW - d'^U^) + ig ) $ - r (9^$)] , (18) 

using Eqs. (2) and (5)~(8). Imposing the variational principle of action and the 
Euler-Lagrange equation in term of $ [42], 

— -a^^ = o (19) 

we find the EOM for a single matter as follow, 

{d^ + ml + 2g^U^)<!> + ^X<^>^ = . (20) 

for a real $ field. 

This result leads to a solitonic wave equation for A 7^ described by the well- 
known nonlinear Klein- Gordon equation, 

{d''+ml)^ + ^X^' = , (21) 

with m| = m| + 2g'^U'^, and C/^ = (fP from Eqs. (9) and (10). Here A determines 
the 'level of non-linearity' for the Klein-Gordon equation. If one puts A ~ m|. we 
arrive at the sine-Gordon equation in 4— dimensional space-time (t, x), — 9^$ — 
m| sin $ = using sin $ ~ $ — -|- • • •. This kind of equation often appears in the 
models based on the coupled pendulum chains pioneered by Englander et.al. [12]. 
However, we should note that the equality A ~ m| in this model doesn't make sense 
since A and m| have different dimensions. In this paper, rather than considering 
that special case, let us solve Eq. (21) in a general way. 

For the sake of simplicity, we consider a traveling wave in 2— dimensional space- 
time {t,x), i.e. ^{x') = ^{x — Ct), where C is a phase velocity. Since d^^ — C^9^,$ 
and = 9^,$, Eq. (21) can be rewritten as, 

a2,$ + m|$ + A$^ = , (22) 

with ml = m|/(C^-l) and A = A/[3!(C^ - 1)]. Assuming that = v is a. 
constant makes m$ to also be a constant. Hence we can multiply both sides of Eq. 
(22) with 9a;/ $ to obtain. 



(5,,$)Vm|$2 + ^A$^ 



. (23) 
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Figure 1: The solitonic wave function for a 2-dimensional DNA as a function of x' 
with the couphng constants (7 = 0.1 (sohd hne) and g — ^ (dashed hne) for a fixed 
parameter set (m$, 0, C, A) = (1, 1, 2, 4). 

Concerning that the quantum wave function $ has the Gaussian distribution, it is 
integrable and then leads to the foUowing differential equation, 



Through standard mathematical procedures, we can straightforwardly get the solu- 
tion. 



for A > 0, or |C| > 1. 

The non-relativistic limit can be obtained by performing a transformation t — > 
T = it in Eq. (21), and putting 7—^1 respectively. This leads to the same result as 
Eq. (25), but t is replaced with —it. The behavior of this solitonic wave function 
is depicted in Fig. 1 as a function of x' with (solid line) and without (dashed 
line) surrounding fluid for a fixed parameter set. Anyway, the fluid contribution is 
independent on its velocity since the effective mass m| is shifted by U"^ — (f?. From 
the figure, we can conclude that the large and small amplitudes can be considered 
as the effects of fiuid surrounding the DNA. 

3.2 Double stranded DNA : the non-Abelian SU(2) model 

Now let us apply the present lagrangian in a more realistic case of double stranded 
DNA. Concerning the smallest group beyond U(l), we take the SU(2) group to 




(24) 




(25) 
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construct the model. In this group, we have 2 sub-matters in a doublet of matter 
field as Eq. (3) with n = 2. 

Since we have only 2 different states of matter, $i and $2, it is convincing to split 
the nucleotide to be a phosphate and a nucleoside consisting of sugar and base. So, 
the interaction between two nucleotides, which further form the backbone of DNA 
molecule, is attributed to the interaction of two different matters, i. e. phosphate and 
nucleoside. On the other hand, the base pair is revealed as the interaction between 
two identical matters, i. e. two neighboring nucleosides belonging to different strands. 
The model is schematically illustrated in Fig. 2 where we have assigned $1 for 
the nucleosides and $2 for the phosphates. Following the allowed interactions in 
Eqs. (11)~(17), we can easily estimate the order of magnitudes for each interaction 
relevant to Fig. 2 as listed in Fig. 3. 

From Fig. 3, it is straightforward to deduce that /i bound is materialized by 
vertex A, while vertex B is responsible for I2 and Is bounds. Anyway, we should note 
that there are another diagrams with spring loops in the vertices A and B, however 
they would be vanishing due to the anti-symmetric /"^'^. Now, we unfortunately 
face a problem on distinguishing I2 with 1^ in Fig. 2. It is quite natural to consider 
I3 must be larger than I2, since the backbone should be rather strongly tied and 
rigid than the nucleotide. Therefore in order to resolve this problem we propose an 
additional contribution to coming from interacting fluid (either fluid absorption or 
emission) with matters depicted in vertex D of Fig. 3. Of course, so Ji could receive 
additional contribution from vertex C too. This scenario could be understood in 
the following way. Since the backbone is more open to surrounding fluid than the 
phosphate— nucleoside encaged in the nucleotide bound-state, its interaction with 
surrounding fluid would contributes more signiflcantly, and then should be taken 
into account. 

We might remark that in the present case the nucleoside, consisting of sugar and 
base, should be considered as a well-conflned bound-state. So we are not going into 
insight to investigate its structure. In consequence of this, we can not distinguish the 
A— T (adenine— thymine) with the G— C (guanine— cytosine) base pairs. Although 
in principle, these might be explained using multi-loops gauge boson exchanges 
inside nucleosides, and two different base pairs are attributed to the ffuid velocities 
in the fluid loops (the second diagram of vertex A in Fig. 3) with opposite sign. 
However we postpone this point in this paper since it would require larger group 
like SU(3) containing more matter fields. Anyway, the opposite torsional motions 
of neighboring strands forming a DNA molecule can be explained, at time being, 
qualitatively by considering the surrounding fluid in both strands have the same 
velocities (v) but with opposite sign each other. 

Now, we investigate the EOM in SU(2) as done in Sec. 3.1. Substituting the full 
lagrangian, Eqs. (5)~(7), into Eq. (19), we obtain for each element of matter, 

a2$i+m|$,-2^ (a^t/|')$2 + ^A ($^ + $2)$,-4^t/A'(a^$2) = 0,(26) 

92^2 + ml $2 + 2g (d^U^) $1 + ($? + $3 $2 + ^9 (d^^i) = , (27) 
for real flelds $i (i : 1,2). Eqs. (4), (26) and (27) immediately yield the following 
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sugar Q phosphate Q base nucleoside 



Figure 2: The double stranded DNA in the non-Abehan SU(2) model with nucle- 
osides and phosphates are represented by $i and $2- The vertices Ji, I2 and I3 
denote different types of interactions connecting nucleosides ($1 — $1) manifest- 
ing base pairs in neighboring strands, nucleoside— phosphate ($1 — $2) within a 
nucleotide, and nucleoside— phosphate ($1 — $2) between nucleotides in a strand. 



EOM, 

(92 + _ a2U^d^) <^> + ^ A ($^$$) = , (28) 

for constant fluid velocity and 4>. Comparing this result with Eq. (21), contribution 
from the interacting fluid medium also appears in the third term but it contributes 
differently. Using Eqs. (9) and (10) we arrive at non-relativistic hmit, 

+ dl^ - m|$ + 4ga2 (p (9^$ T « v • d^^) - ($$^*) = , (29) 

for V2 = V. This is the nonlinear EOM governing the double stranded DNA dynam- 
ics with surrounding ffuid medium in the present theory. The plus and minus signs 
show the dynamics of a strand and its counterpart surrounded by the ffuids with 
opposite velocities. 

Obviously, in contrast with the U(l) case it is hard to solve Eq. (29) exactly. 
For the sake of simpliflcation, let us consider 2-dimensional {t,x) case of Eq. (28), 

where at = —4:iga2j(f) and ax — Aiga2jVx(f)- Borrowing the traveling wave ^{x') = 
^{x — Ct) as before we obtain, 

|£-a§ + *|4.A*3^0. (31) 
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0(g) 



0(g) 



Figure 3: The 2— point interactions and its first order contents relevant for double 
stranded DNA in Fig. 2 and its order of magnitudes. The plain and spring lines 
indicate matter and fluid fields. 



with a = {Cat + ctx)/{C^ - 1), ml = m|/(C2 _ ^ = a/(3!(C2 - 1)). For 

A = it coincides with the equation of inharmonic oscillator, Eq. (25). 

For further simplification, wc assume that A is small enough such that the last 
term in Eq. (31) can be treated pcrturbativcly. Then, we can expand the mass m$ 
in term of A, i.e. rh^ ~ m$,) + Arn$^, and $ ~ $o + A$i up to 0(A) accuracy. Now 
we are ready to solve Eq. (31) order by order. 

The lowest order with respect to A, i.e. 0(1), satisfies the following equation. 



(32) 



dx'^ dx' 

Following the standard mathematical procedures the solution can in general be 



•^•'^"^^^^ ^^^^^^^^^^^^^^^^^ ^^^^^^^^^^ ^^^^ ^^^^^^^^^ 

expressed in the form of $o — N^e^+^' + N^c'^-^' with k± = ^ 



Therefore, the solution of Eq. (32) is simply, 

1 

cosh 



for 



$0 = 2A^oe^^''x < 



sm 




1 - ^^0^ 



ml, - ^x' 



for > 4m||j 



for < 4m|^ 



(33) 



by putting the normalization factor to be A^"q 



Nq. Each solution is corre- 



sponding to the over-damped, damped and regular harmonic oscillators respectively. 
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Subsequently, the next order, i.e. 0(A), is governed by the equation, 

^^-tt^ + m|^*i + 2m$>$,$o + $o = 0- (34) 
The over-damped $o in Eq. (33) yields the general solution for $i should be, 

$1 = AT+et^-' + iVfe^^' . (35) 



Substituting the over-damped $o S'-iid Eq. (35) into Eq. (34), we obtain A^^ = 
-(32iVo3)/(3a2 + 4m|J. In the present case the first term in Eq. (35) is vanishing 
for any since = 4m|jj. This also leads to the result = since mq,^, Nq 0. 
Finally, 

32N^ 



$1 = 



3q;2 -I- 4m|^ 



3a I 

e 2 



(36) 



and m$ = m$,). 

In the second case of damped $o, we make use of the equality cosh^ x = 1/2 cosh(3a;)-|- 
coshx to obtain. 



dx' 



52 



^2 ^/ 



+4iVn"e 



^2 ^"')cosh| ^/^-ml^a;' 



. (37) 



Since coshx < cosh(3a;) and the last term is enhanced only by a factor of as small 
as 2, Eq. (37) can be approximately reduced to be. 



(38) 



The solution is given by. 



$1 = e 2 



7Vi+ cosh 3 



- - mi^x I + TVi sinh I 3 y — - 



— mi a;' 



Again substituting it into Eq. (38) yields, 

9 



—N; 



- 3m|^j 



4" ^ 3^4 + 36m|^ - 20a'^ml^ ' 
9 



a\ oP' — Ami 



$0 



(39) 



(40) 



(41) 



From these results, for a;' > we can safely omit the sub-dominant sinh term in Eq. 
(39), also because > since oc^ > 4m|^. Hence, 



9^0 - 3m| 

$1 = — — 

4 ° 3a4 + 36m|^ - 20a2rfi|^ 



e^a"^' cosh f 3^/^^ rn^^^^x 



^2 I 



(42) 
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Figure 4: The wave function for 2-dimensional double stranded DNA as functions 
of x' with = 1.9 (sohd hue) and v.j. = 1.7 (dashed hne) for a fixed parameter sets 
(m$„,0,C,A,^,iVo) = (l,l,2,4,l,l). 



The last case of regular harmonic oscillator is governed by the following equation, 



a- 



dx' 



1 , j:^2 



m|^$i - 2A^o^e^^ sin ( 3\/m|„ - —x' 



+2A^o^e^^ 3 + 



sm 



(52 



^~4o--r^' = 0' (43) 



using the relation 4sin^a; = Ssinx — sin(3a;). In contrary with previous cases the 
general solution for Eq. (43) is complicated. So, let us assume here that the more 
rapid oscillation term, i.e. the fourth term, is dominant than the last one which 
reduces the equation to be. 



- a—j + m|^$i - 2N^e 2 ^ sin ?>\ - 



dx'^ dx 
Hence the general solution is simply. 



ny2 



a 



-X 







$1 



3a / 
e 2 ^ 



q;2 



cos 3\ - —x' + sin 3\ ml^ - —x' 



(44) 



(45) 



Following similar procedures as before, 

4 



TV- = 



a 4 /4ml — 

IN^ ^ - 

3 ° 3q;4 + 36m|^ - 1652^1^ ' 

3 ° 354 + 36m|^ - 16a^ml„ ' 



(46) 
(47) 



^0 
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In non-relativistic case, by definition the condition m$o > \2g{C — Vx)(l)\ should be 
fulfilled. Obviously, for large enough (a^ <C 4:rh%^ or m^^ 3> \2g{C — Vx)4>\) the 
solution is dominated by term, while both terms are comparable for 4m|jj 
or m$o — >■ \2g{C — Vx)(l)\. These arguments lead to the result. 




We should remark that up to the current accuracy there is no need in all cases to 
calculate the leading order of mass, . As a typical example, the wave function 

for the harmonically oscillating, i.e. the sum of Eqs. (33) and (48), double stranded 
DNA is depicted in Fig. 4 as a function of x' for certain velocities. It can also be 
seen that the oscillation is sensitive to the fiuid velocity. 



4 Summary and discussion 

We have introduced a new type of model to describe DNA using the gauge invariant 
fiuid lagrangian. The lagrangian is able to accomodate various internal motions of 
DNA, from the single bulk to the double stranded of DNA as done in the preceeding 
section. The EOM's and its solutions for two typical cases using the Abelian U(l) 
and non-Abelian SU(2) lagrangians have been derived and investigated. 

In the case of Abelian U(l) lagrangian, we have seen from Eq. (25) that the 
interacting fiuid medium characterized by the coupling constant g infiuences the 
magnitude and also the width (associated to the dispersion or steppening rate) of 
solitonic wave equation as well, but regardless with the fiuid velocity. On the other 
hand, obviously the matter self-interaction represented by its coupling constant A 
could change only the magnitude and not the dispersion or steppening rate of soliton. 
Actually, this provides a natural explanation for small and large amplitude regions 
of the internal motion of a single bulk of DNA immersed in bio-fiuid without adding 
any new terms by hand as done in some previous works [43]. Furthermore, that 
contribution shifts the matter mass ni<j, to be m$. This is the so-called running 
mass induced by the dynamical fluctuation of internal kinematics in the system 
as a result of interaction between matter and fluid. However, the result is again 
independent on the fiuid velocity. 

In the second case, using the non-Abelian SU(2) lagrangian we have constructed a 
model for double stranded DNA in detail up to the level of its constituents, except for 
sugar and base composing the nucleoside. It has been shown that the EOM follows 
a similar form as in the U(l) model, but the interacting fluid medium contributes 
in different way. The model requires that the DNA polymer would exist if and only 
if it resides in a fiuid medium, represented by I2 and I3 bounds realized by fiuid- 
matter interactions. Otherwise, the binding interactions I2 and I3 would vanish and 
the strands are broken. These results could be used to explain the deformation of 
DNA molecules associated with vanishing interactions in /i,2,3- In contrast with the 
previous case, the fiuid velocity plays an important role and changes the dynamics 
drastically, namely the highly damped, damped and regular harmonic oscillators. 
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This supports a conclusion obtained in [30], that is the effect of hydrodynamic 
interactions on the dynamics of DNA translocation depends on the fluid velocity. 

As mentioned earlier, both strands in a double stranded molecule are considered 
to follow the same EOM as Eq. (29) with opposite sign of fluid velocities. In 
contrary, the single fragments of DNA belong to those strands arc governed by Eq. 
(21) independent on the fluid velocity, and should behave identically no matter with 
the directions of its surrounding fluid velocities. 

Further studies can be done using the lattice gauge simulation to calculate 
numerically, for instance the finite temperature partition function density Z — 
exp(l/T J ^xC). This kind of macroscopic ensemble provides direct relation be- 
tween the internal dynamics of DNA and some physical observables like tempera- 
ture and so on. Actually this is the main advantage of deploying the gauge invariant 
lagrangian like the present one. Such numerical calculations would be able to sim- 
ulate quantitatively some phenomena in DNA like critical temperature or pressure 
related to the deformation of DNA molecules, etc. For example, one can investigate 
the critical temperature as a double stranded DNA is splitted into single strands 
[44], i.e. /i ^ in the present model. Such studies are in the progress. 
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